arXiv:1509.03574v2 [cs.DM] 18 Sep 2015 


On extremal trees with respect to the f’-index 


20 Juli 2015 



°‘Institut fiir Informatik, Freie Universitdt Berlin, 
Takustrafie 9, D-14195 Berlin, Germany 

E-mail: abdo@mi.fu-berlin.de 

^Hochschule fiir Technik und Wirtschaft Berlin, 
Wilhelminenhofstrafie 75A, D-12459 Berlin, Germany 

E-mail; darko. dimitrovl l@gmail. com 

'^Faculty of Science, University of Kragujevac, Kragujevac, Serbia 

E-mail: gutman@kg. ac. rs 

‘’'State University of Novi Pazar, Novi Pazar, Serbia 


Abstract 


In a study on the structure-dependency of the total vr-electron energy from 1972, 
Trinajstic and one of the present authors have shown that it depends on the sums 
Yhv&v cind Yhv&v where d{v) is the degree of a vertex v of the underling 

molecular graph G. The hrst sum was later named first Zagreb index and over the 
years became one of the most investigated graph-based molecular structure descrip¬ 
tors. On the other hand, the second sum, except in very few works on the general 
hrst Zagreb index and the zeroth-order general Randic index, has been almost com¬ 
pletely neglected. Recently, this second sum was named forgotten index, or shortly 
the F-index, and shown to have an exceptional applicative potential. In this paper 
we examine the trees extremal with respect to the F-index. 


^Corresponding author 



1 Introduction 


The first and the second Zagreb indices, introduced in 1972 16 , are among the oldest 


graph-based molecular structure descriptors, so-called topological indices ll8,ll 


For a graph G with a vertex set V{G) and an edge set E{G), these are dehned as 
Mi(G')= d{vf and M 2 {G) = E d{u)d{v), 

v£V(G) uv£E{G) 

where the degree of a vertex v G V{G), denoted by d{y), is the number of the hrst 
neighbors of v. 

Over the years these indices have been thoroughly examined and used to study 
molecular complexity, chirality, ZE-isomerism and hetero-systems. More about their 
physico-chemical applications and mathematical properties can be found in [Io|[ITp4] 
and |6| [I^[l4pH[^ , respectively, as well as in the references cited therein. The sum of 
squares of vertex degrees was also independently studied in quite a few mathematical 
papers 


25 


In an early work on the structure-dependency of the total vr-electron energy 16 


beside the first Zagreb index, it was indicated that another term on which this energy 
depends is of the form 

F(G)= d(v)K 

v&V{G) 

For unexplainable reasons, the above sum, except (implicitly) in a few works about 


the general first Zagreb index 19,20 and the zeroth-order general Randic index 17 


has been completely neglected. Very recently, Furtula and one of the present authors 
succeeded to demonstrate that F{G) has a very promising applicative potential [^. 
They proposed that F{G) be named the forgotten topological index, or shortly the 
F-index. 

Before the publication of the paper [^, the F-index was not examined as such. On 
the other hand, in some earlier studies on degree-based graph invariants, it appears 
as a special case. 

The general first Zagreb index of a graph G is defined as 
Mf{G)= ^ ^ [d{u)°‘~^ + d{v)^~^], for a G M, a 7 ^ 0, a 7 ^ 1, 

v£V{G) uveE{G) 












and its first occurrence in the literature seems to be the work 19 by Li and Zhao from 


2004. Observe that (G) = F{G). In 19 the trees with the hrst three smallest and 


largest general hrst Zagreb index were characterized. In 19,20 , among other things, 


it was shown that for a > 1, the star S'„ is the tree on n vertices with maximal M" 
whereas the path is the tree on n vertices with minimal M“-value. Needless to say 
that these results directly apply to the P-index (a = 3). More results and information 
about the general hrst Zagreb index can be found in . 

The Randic (or connectivity) index was introduced by Randic in 1975 and is 
dehned as 

R = R(G) = y , ^ . 

uv^G) Vd{u)d{v) 

Later in 1977, Kier and Hall have introduced the so-called zeroth-order Randic 
index 

°P = °P(G)= y d(n)"5, 

v&V{G) 

which was generalized in 2005 by Li and Zheng [^. It was named the zeroth-order 
general Randic index and was dehned as 

°p„ = op„(G)= y d{vr 

v&V{G) 

for any real number a. Evidently, ^R 2 ,{G) = F{G). 


In 17 graphs of maximal degree at most 4 (molecular graphs ), with given number 


of vertices and edges, and with extremal (maximum or minimum) zeroth-order general 
Randic index were characterized. This, again, in the special case a = 3 renders results 
for the P-index. 

Here, we extend the work on trees with extremal values of P-index, by considering 
trees with bounded maximal degree. Since the paths are the trees with minimal P- 
index also in this case, we consider here the characterization of trees with bounded 
maximal degree that have maximal P-value. 

In the sequel we introduce notation that will be used in the rest of the paper. 
For u,v,x,y G V{G) such that uv G E{G),xy ^ E{G), we denote hy G — uv the 
graph that is obtained by deleting the edge uv from G and hy G + xy the graph 
that is obtained by adding the edge xy to G. By A = A(G) we denote the maximal 










degree of G. A sequence D = [di, ^ 2 , • • •, is graphical if there is a graph whose 
vertex degrees are d*, i = 1,..., n. If in addition di > d 2 > ■ ■ ■ > dn, then D is 
a degree seguence. Furthermore, the notation D{G) = , ■ ■ ■ means that 

the degree sequence is comprised of n* vertices of degree Xi, where i = 1 , 2 , • • ■ ,t. 

A tree is said to be rooted if one of its vertices has been designated as the root. 
In a rooted tree, the parent of a vertex is the vertex adjacent to it on the path to 
the root; every vertex except the root has a unique parent. A vertex is a parent of a 
subtree, if this subtree is attached to the vertex. A child of a vertex u is a vertex of 
which V is the parent. 

2 Results 

First, we characterize trees with maximal degree at most A that have maximal F- 
index. 

2.1 Trees with bounded maximal degree 

Theorem 2.1. Let T be a tree with maximal F-index among the trees with n vertices 
and maximal degree at most A. Then the following holds: 

(i) // (n — 2) mod (A — 1) = 0, then T contains vertices of degree A and 
vertices of degree 1 . 

(a) Otherwise, T contains vertices of degree A, vertices of degree 

1 and one vertex of degree x, where x is uniquely determined by 2<x<A — 1 
and {n — 1 — x) mod (A — 1 ) = 0. 

Proof. We may assume that T is a rooted tree, whose root is a vertex with degree A. 
First, we show that T is comprised of vertices of degrees A and 1, and in some cases 
of one additional vertex u, with 2 < d{u) < A — 1 . Assume that this is not true and 
that T has more than one such vertex. Denote by lA the set of all vertices of T whose 
degrees are different from A and 1. We assume that for the degrees of the vertices of 
Vd the order d{u) > ■ ■ ■ > d{v) holds. Let tc be a child vertex of v. 



Delete the edge vw and add the edge uw to T, obtaining a tree T'. After this 
transformation, the degree of u increases by one, while the degree of v decreases by 
one, and the degree set Vd changed into V^. It holds that 

F(T') - F{T) = {d{u) + 1)^ - d{uf + {d{w) - 1)^ - d{wf > 0. 


It holds that |V^'| < |Vd|, with strict inequality if d{u) = A —1 or d{v) = 2. If |V^'| > 2 
we choose from I/j a vertex with maximal and a vertex with minimal degree, and 
we repeat the above operation, obtaining a tree with larger F-index. We proceed 
on iteratively with the same type of transformation until the transformed Vd has 
cardinality 1 or 0, and thus obtain a contradiction to the assumption that T has 
more that one vertex with degree different from 1 and A. It follows that T has a 
degree sequence or 1”^], 2 < a; < A — 1. 

Observe that any transformation on T will result in a tree with \Vd\ > 1 or will 
disconnect T. Thus, we conclude that T must have one of the two above presented 
degree sequences. 

Next, with respect to the cardinality of |Vd|, we determine the parameters tta , 
rii , and x. 


Case 1. \Vd\ = 0. In this case the degree sequence of T is [A1*^1] where the 
equations 


A A 

Ui = + ni = n and i ni = An^ + ni = 2(n 

i=\ i=l 

hold. The above equations have the integer solution 


Ua 

for (n — 2) mod (A — 1) 


n — 2 
A- 1’ 

0 . 


Til = 


n(A-2) + 2 
A - 1 


1 ) 


Case 2. \Vd\ = 1. Here the degree sequence of T is [A"'^,a;^, with 


n* = ua + ni + 1 = n 


and 


i rii = Aua + ni + a; = 2(n — 1). 


i=l 


2=1 


The above equations give the integer solution 


ua = 


n — 1 — X 
A - 1 ’ 


ni = 


(n — 1)(A — 2) + a; 
A- 1 ’ 


for (n — 1 — a;) mod (A — 1) = 0. 


□ 



As straightforward consequence of Theorem 2 . 1 , we obtain the maximal value of 
the F-index for trees with maximal degree A. 

Corollary 2.1. Let T be a tree with maximal F-index among the trees with n vertices 
and maximal degree at most A. Then, 

(i) if {n — 2) mod (A — 1) = 0, 


F{T) = A(A + l)(n - 2) + 2(n - 1); 


(a) otherwise, 


F{T) = (A^ + A + 2)(n - 1) - (A^ + A + l)x + x^, 


where X is uniguely determined by2 < x < A —1 andn — l—x mod (A —1) = 0. 


Proof. If (n —2) mod (A —1) = 0, by Theorem 2.1, T has degree sequence [A”^, 1"^], 
where ua = (n — 2 )/(A — 1 ) and rii = (? 7 ,(A — 2 ) + 2 )/(A — 1 ). Thus, 


F{T) = A 


3 n — 2 ^ ?7,(A — 2) + 2 


+ 


= A(A + l)(n-2) + 2(n-l). 


A - 1 A - 1 

Otherwise, T has degree sequence [A"'^,x^, 1"^], where ua = (n — 1 — x)/{/S. — 1), 
Hi = ((n — 1)(A — 2) + a;)/(A — 1). Then, 


F{T) = A 


3 n — (a; + 1) (n — 1)(A — 2) + a; 3 


A - 1 


+ 


A - 1 


+ a;" 


(A3 + A - 2)n + (1 - A3)a; - (A^ + A - 2) 3 

“ A^I 

= (A^ + A + 2)(n — 1) — (A^ + A + l)a; + x^. 


□ 


2.2 Molecular trees 

For the special case of molecular trees, i.e., A < 4, we obtain the following result. 
Recall that such trees provide the graph representation of the so-called saturated 
hydrocarbons or alkanes (T, 15 and are of major importance in theoretical chemistry. 


Theorem 2.2. Let T be a molecular tree with maximal F-index among the trees with 
n vertices. Then the following holds: 











If {n —2) mod 3 = 0, then T contains vertices of degree 4 and vertices 
of degree 1. Its F-index is 


F{T) = 22n - 42. 

Otherwise, T contains vertices of degree A, vertices of degree 1 

and one vertex of degree x, where x is uniguely determined by 2 < x < 3 and 
{n — 1 — x) mod 3 = 0. Its F-index is 

F{T) = 22(n- 1) -2lx + x^. 


Proof. For n = 2,3,4, 5 the star Sn maximizes the F-index . For n = 6 and A < 4, 
the possible degree sequences of T are (4,2,1,1,1,1), (3, 3,1,1,1,1), (3, 2,2,1,1,1) 
and (2, 2, 2,2,1,1). The hrst of these degree sequences corresponds to the largest 
F-index, F(T). 

Thus, the corresponding degree sequences of T for n = 2,3,4, 5,6 satisfy the 
theorem. 

In the rest of the proof we assume that n >7. First, we show that T has maximal 
degree A = 4. A tree with maximal degree A we denote by T^,. By Corollary 2.1[ z), 
F{T 2 ) = 8n — 14. Further, we distinguish between two case. 

Case 1. {n — 2) mod 2 = 0. 

By Corollary 113 i), we have that F{Ts) = 14n — 26. Since (n — 2 ) mod 2 = 0, it 
follows that {n—2) mod 3 7 ^ 0, and by Corollary |2.1 [ m), F(T 4 ) = 22(n —1) —21x+x^, 
where x = 2 or a: = 3. For n>7, F(T 4 ) > F{Ts) > F(T 2 ) is satished. 

Case 2. {n — 2) mod 2 7 ^ 0. 

By Corollary 2.l\ ii), F{Ts) = 14(n — 1) — 13x + x^, where x = 2 or x = 3. 

Subcase 2.1. (n — 2) mod 3 = 0. 

By Corollary |2.1[ i), FifT.^) = 22n — 42. For n > 7 and x = 2,3, it holds that 
F(T 4 ) > F{T,) > F{T,). 

Subcase 2.2. (n — 2) mod 3 7 ^ 0. 

By Corollary 2.l\ ii), FifF/f) = 22 (? 7 ,—1)—21x+x^. Again, a straightforward calculation 
yields that F(T 4 ) > F{T^) > F(T 2 ). 
















So, we have shown that for n > 3,4, 5 the tree with maximal F-index has maximal 
degree n — 1, and for n > 5 it has maximal degree A = 4. Setting A = 4 in 


Theorem 2.1 and Corollary 2.1, we complete the proof. 


□ 


We would like to note that the results of Theorem 12.21 coincide with the results of 


Theorem 2.2 in 17 , pertaining to the zeroth-order general Randic index of a graph, 
when a = 3 and m = n — 1. 

In the sequel we present some computational results obtain by integer linear pro¬ 
gramming and by exhaustive computer search of trees with maximal F-index. 


2.3 Computational results 

Degree sequences of trees with bounded maximal degree and maximal F(T) can be 
completely described by solving the integer linear programming problem in Lemma 
2.1[ We denote the number of vertices of tree T of degree i by n^, Vi = 1, 2, • • • , A 
and the number of edges, connecting a vertex of degree i with a vertex of degree j, 
by rriij. 

Lemma 2.1. Let T be a tree with maximal F-index among all trees with n vertices 
and maximal degree A. Then, the trees with maximal F-index are completely described 
by solving the maximization problem 

A 

maximize F = i^ ui 

i=l 

A 

subject to rij = n, 

2=1 

A 

'^iui = 2{n- 1 ), 

2=1 

A-1 

^ Ui < 1, 

i=2 
A 

mij -|- 2mjj = iui 1 < i < A, 

0<?7,j<n — 1, l<i<A, 

0 < niij < n — 2, 


1 < *, J < A. 









For experimental proposes we first consider the molecular trees of orders n = 
4, • • • , 20. In this case, the previous integer linear programming model can be repre¬ 
sented as follows. 


s. t. 


F = Til 

+ 

00 

to 

+ 

27n3 

+ 

64^4 



ni 

+ 

n2 

+ 

ns 

+ 

77-4 


= n, 

Til 

+ 

2n2 

+ 

3^3 

+ 

4?7,4 


= 2{n 



n2 

+ 

ns 




^ 1, 

-ni 

+ 

mi2 

+ 

mis 

+ 

mi4 


= 0, 

-2n2 

+ 

mi2 

+ 

2m22 

+ 

m23 

+m24 

= 0, 

-3n3 

+ 

mi3 

+ 

m2s 

+ 

2mss 

+m34 

= 0, 

—4n4 

+ 

mi4 

+ 

17124 

+ 

77134 

+2m44 

= 0, 

0 ^ Uj ^ 

n — 

o 

/A 

/A 

-2, 

1 < i 

1 

VI 

VI 

VI 


( 1 ) 


The above integer linear program was implemented in Matlab 23 , and run on 2.3 


GHz Intel Core i5 processor with 4GB 1333 MHz DDR3 RAM. Beside the degree 
sequence that maximize the F-index, it gives only one corresponding tree. However, 
for a given degree sequence there usually exist several non-isomorphic trees with the 
same degree sequence. In view of this, the degree sequence as well as all corresponding 
trees were obtained by exhaustive computer search using the mathematical software 


Sage 26 . In the contrast to the integer linear program that gives a solution for n of 


order 1000000 and A < 400 in less than a minute, the exhaustive computer search 
with Sage for n = 20 took several hours. On the other hand, due to the closed-form 


solutions in Theorem |2.1[ we obtain the degree sequence that maximizes the F-index 
for arbitrary n and arbitrary A in few milliseconds. 

The experimental results for A = 4 and A = 5 are given in Tables [T] and 
respectively, while their corresponding trees are given in Figures [T] and 







n 

Sage 

LP solution - Matlab 

p 

D(T) 

#T 



Non-zero 

variables 


± 

4 

[3M3] 

1 

Hi 

3 

1 


"^1,3 

3 



30 

5 

[4^,1^] 

1 

ni 

4 

n4 

1 


7771,4 

4 



68 

6 

[ 41 , 21 , 1 ^] 

1 

ni 

4 

n2 

1 

rii 

1 

7771,2 

1 

777i,4 

3 

7772,4 

1 

76 

7 

[ 41 , 31 , 1 ^] 

1 

ni 

5 

1 

774 

1 

"tl,3 

2 

7771,4 

3 

?^3,4 

1 

96 

8 

[ 42 , 16 ] 

1 

ni 

6 

Hi 

2 


777i,4 

6 

7774,4 

1 ’ 


134 

9 

[ 42 , 21 , 16 ] 

2 

ni 

6 

n2 

1 

774 

2 

7771,4 

6 

7772,4 

2 


142 

10 

[ 42 , 31 , 1 ^ 

2 

ni 

7 

1 

774 

2 

^^1,3 

1 

777i,4 

6 

"^3,4 

2 

162 

11 

[ 43 , 18 ] 

1 

ni 

8 

n4 

3 


7771,4 

8 

7774,4 

2 ’ 


200 

12 

[ 48 , 21 , 18 ] 

3 

Hi 

8 

n2 

1 

774 

3 

777i,4 

8 

7772,4 

2 

7774,4 

1 ’ 

208 

13 

[ 48 , 31 , 11 *] 

4 

ni 

9 

nz 

1 

774 

3 

7771,4 

9 

3 


228 

14 

[4^ 116] 

2 

ni 

10 

714 

4 


777i,4 

10 

7774,4 

3 


266 

15 

[4^2l,ll0] 

6 

ni 

10 

n2 

1 

774 

4 

7771,4 

10 

7772,4 

2 

7774,4 

2 

274 

16 

[4^3l,lll] 

8 

ni 

11 

1 

774 

4 

777i,4 

11 

mzA 

3 

7774,4 

1 ’ 

294 

17 

[45, 112] 

3 

ni 

12 

774 

5 


7771,4 

12 

7774,4 

4 ’ 


332 

18 

[ 46 , 21 , 112 ] 

14 

ni 

12 

n2 

1 

774 

5 

777i,4 

12 

7772,4 

2 

7774,4 

3 ’ 

340 

19 

[ 46 , 31 , 118 ] 

17 

ni 

13 

nz 

1 

774 

5 

7771,4 

13 

^^3,4 

3 

7774,4 

2 

360 

20 

[ 46 , 114] 

5 

ni 

14 

774 

6 


7771,4 

14 

7774,4 

5 ’ 


398 


Table 1: Extremal molecular trees of order up to 20. First column: the order of a tree. Second 
column: the degree sequences that maximize the E-index and the number of corresponding trees, 
obtained by Sage. Third column: solutions of an integer linear programming problem with one 
realization of each degree sequence, obtained by Matlab. Fourt column: the value of the E-index. 























n 

Sage 



LP solution - 

Matlab 



E 

D(T) 

#T 



Non- 

■zero variables 




4 

[3M^] 

1 

ni 

3 

nz 

1 


"^1,3 

3 




30 

5 

[4^l1 

1 

ni 

4 

rii 

1 


mi,4 

4 




68 

6 


1 

ni 

5 

ns 

1 


^^1,5 

5 




130 

7 

[5',2i,l5] 

1 

ni 

5 

n2 

1 

ns 

1 

mi,2 

1 

^^1,5 

4 

tn2,5 

1 


138 

8 

[S', 31 , 16 ] 

1 

ni 

6 

nz 

1 

ns 

1 

"^1,3 

2 

^^1,5 

4 

tn3,5 

1 


158 

9 

[ 51 , 41 , 1 ^ 

1 

ni 

7 

rii 

1 

ns 

1 

mi,4 

3 

4 

tn4,s 

1 


196 

10 

[5^ 1«] 

1 

ni 

8 

ns 

2 


^^1,5 

8 

^^5,5 

1 



258 

11 

[ 52 , 21 , 18 ] 

1 

rii 

8 

^2 

1 

ns 

2 

^^1,5 

8 

m 2 ,5 

2 



266 

12 

[52,3i,li»] 

2 

ni 

9 

riz 

1 

ns 

2 

"^1,3 

2 

^^1,5 

7 

tn3,5 

1 

tns,5 

1 

286 

13 

[ 52 , 41 , 110 ] 

2 

ni 

10 

rii 

1 

ns 

2 

mi,4 

2 

?^1,5 

8 

tn4,s 

2 


324 

14 

[ 58 , 111 ] 

1 

ni 

11 

ns 

3 


^^1,5 

11 

^^5,5 

2 



326 

15 

[58,21,111] 

3 

ni 

11 

^2 

1 

ns 

3 

^^1,5 

11 

m 2,5 

2 

tns,5 

1 


394 

16 

[ 58 , 31 , 112 ] 

4 

ni 

12 

riz 

1 

ns 

3 

"^1,3 

2 

"^1,5 

10 

tn3,5 

1 

tns,5 

2 

414 

17 

[ 58 , 41 , 118 ] 

4 

ni 

13 

rii 

1 

ns 

3 

mi,4 

1 

trips 

12 

tn4,5 

3 


452 

18 

[5^ 11^] 

2 

ni 

14 

ns 

4 


^^1,5 

14 

ms ,5 

3 



514 

19 

[5^2l,ll^] 

7 

ni 

14 

^2 

1 

ns 

4 

mi,2 

1 

trips 

13 

tn2,5 

1 

tns,5 

3 

522 

20 

[5^3l,ll8] 

8 

ni 

15 

1 

ns 

4 

^^1,3 

2 

trips 

13 

tn3,5 

1 

tns,5 

3 

542 


Table 2: Extremal trees of order up to 20 and maximal degree A = 5. First column: the order 
of a tree. Second column: the degree sequences that maximize the E-index and the number of 
corresponding trees, obtained by Sage. Third column: solutions of an integer linear programming 
problem with one realization of each degree sequence, obtained by Matlab. Fourt column: the value 
of the E-index. 





n = 20 


Figure 2: Extremal trees of order n = 4, 5, • • • ,20 and A = 5 with respect to the F-index. 
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